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Abstract 

The isotope effect in the scattering of methane is studied by wavepacket simu- 
lations of oriented CH4 and CD4 molecules from a flat surface including all nine 
internal vibrations. At a translational energy up to 96 kJ/mol we find that the 
scattering is still predominantly elastic, but less so for CD4. Energy distribution 
analysis of the kinetic energy per mode and the potential energy surface terms, 
when the molecule hits the surface, are used in combination with vibrational exci- 
tations and the corresponding deformation. They indicate that the orientation with 
three bonds pointing towards the surface is mostly responsible for the isotope effect 
in the methane dissociation. 



keywords: Computer simulations, Models of surface chemical reaction, Alka- 
nes, Low index single crystal surfaces 



1 Introduction 

The dissociation of methane on transition metals is an important reaction in 
catalysis. It is the rate limiting step in steam reforming to produce syngas. [1] 
It is also prototypical for C-H activation in other processes. A large number of 
molecular beam experiments in which the dissociation energy was measured as 
a function of translational energy have already been done on this system. [2- 
22] These experiments have contributed much to our understanding of the 
mechanism of the dissociation. Some of them observed that vibrationally hot 
CH 4 dissociates more readily than cold CH 4 , with the energy in the inter- 
nal vibrations being about as effective as the translational energy in inducing 
dissociation. [2-8] A more detailed assessment of the importance of the inter- 
nal vibrations could not be made, because of the large number of internal 
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vibrations. A recent molecular beam experiment with laser excitation of the 
v 3 mode succeeded in measuring a dramatical enhancement of the dissociation 
on a Ni(100) surface, but it is still much too low to account for the vibrational 
activation observed in previous studies and indicates that other vibrationally 
excited modes contribute significantly to the reactivity of thermal samples. [22] 

Wavepacket simulations are being used more and more to study the dynamics 
of this kind of molecule surface reactions. The published wavepacket simula- 
tions on the methane dissociation reaction on transition metals have treated 
the methane molecule always as a diatomic up to now. [23-27] Besides the C-H 
bond and molecule surface distance, a combination of other coordinates were 
included, like (multiple) rotations and some lattice motion. None of them 
have looked at the role of the internal vibrations. Various theoretical stud- 
ies have obtained reaction pathways and barriers for dissociation by DFT 
calculations, [28-37] but they cannot explain the role of the vibrational modes 
in the reaction dynamics either. 

A nice way to study reaction dynamics is the use of isotopes. The most re- 
cent wavepacket simulation on the dissociation probability of CH 4 and CD 4 
showed a semiquantitative agreement with the molecular beam experiments of 
Ref. [5], except for the isotope effect and the extracted vibrational efficacy [27] 
The molecular beam study with laser excitation of the v$ asymmetrical stretch 
mode shows that the incorrect vibrational efficacy is caused by the assump- 
tions in the fit procedure that both stretch modes behaves identical. [22] One of 
the possible explanation of the incorrect isotope effect can be the role played 
by the non-included intramolecular vibrations. 

In a previous paper we reported on wavepacket simulations that we have done 
to determine which and to what extent internal vibrations are important for 
the dissociation of CH 4 .[38] We were not able yet to simulate the dissocia- 
tion including all internal vibrations. Instead we simulated the scattering of 
methane, for which all internal vibrations can be included, and used the results 
to deduce consequences for the dissociation. We used model potential energy 
surfaces (PESs) that have been developed with Ni(lll) in mind, but our re- 
sults should hold for other surfaces as well. At a translational energy up to 96 
kJ/mol we found that the scattering is almost completely elastic. Vibrational 
excitations when the molecule hits the surface and the corresponding deforma- 
tion depend on generic features of the potential energy surface. In particular, 
our simulations indicate that for methane to dissociate the interaction of the 
molecule with the surface should lead to an elongated equilibrium C-H bond 
length close to the surface. 

We have been using the multiconfigurational time-dependent Hartree (MCTDH) 
method for our wavepacket simulation, because it can deal with a large num- 
ber of degrees of freedom and with large grids. [39,40] This method has been 
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applied successfully to gas phase reactions and reactions at surfaces. [41-60] 

In this paper we report wavepacket simulations of CD 4 scattering including 
all internal vibrations for fixed orientations, performed on the same model 
PESs as in our previous paper. [38] Translational motion parallel to the sur- 
face and all rotational motion was neglected. No degrees of freedom of the 
surface were included. Experiments show that coupling with these degrees of 
freedom is dependent on the metal surface. For example, the observed surface 
temperature effect are small on Ni[5], but quite large on Pt[7]. As we are only 
interested in the role of internal vibrations, we have not included degrees of 
freedom of the surface to keep the simulations as simple as possible. We will 
discuss the vibrational excitation and the deformation of the CD4 molecule 
when it hits the surface and compare it with CH 4 . Later on we will look at the 
energy distribution of the kinetic energy per mode and the potential energy 
in some terms of the PES with the elongated equilibrium bond length close 
to the surface for both isotopes. The transfer of translational kinetic energy 
towards vibrational kinetic energy gives an indication about the dissociation 
probability, since vibrational kinetic energy helps in overcoming the dissoci- 
ation barrier. It gives a better idea too about which modes are essential to 
include in a more accurate wavepacket simulation of methane dissociation. 
After that we will discuss the implications of this for the dissociation and give 
a summary with some general conclusions. 



2 Computational details 

2. 1 The Potential Energy Surfaces 

We used for the scattering of CD 4 the same model PESs as we did for CH 4 . 
Since we expressed the PES in mass-weighted coordinated the parameters in 
the PESs for CD 4 differs from CH 4 . We will now give an overview of our 
model PESs and the corresponding parameters for CD 4 . The parameters of 
CH 4 for these PESs were already given in Ref. [38], where also more detailed 
information about our assumptions and contour plots of some cross-section of 
the model PESs can be found. 

The PESs we used can all be written as 

Kotal = Vintra + V Bm f, (1) 

where V^ ntr a is the intramolecular PES and V^ ur f is the repulsive interaction 
with the surface. For the V^ ntra we looked at four different types of PESs. The 
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Vi n tra include also for two types changes in the intramolecular potential due to 
interactions with the surface. 



2.1.1 A harmonic potential 

The first one is completely harmonic. We have used normal mode coordinates 
for the internal vibrations, because these are coupled only very weakly. In the 
harmonic approximation this coupling is even absent so that we can write 

^intra as 

1 10 

Vintra = Harm — X kjXj , (2) 
1 i=2 

the summation is over the internal vibrations, X^s are mass- weighted displace- 
ment coordinates and ki are mass-weighted force constants, (see Table 1 for 
definitions and values); (X 1 is the mass- weighted overall translation along the 
surface normal). [61] The force constants have been obtained by fitting them 
on the experimental vibrational frequencies of CH 4 and CD 4 . [62,63] 

We have assumed that the repulsive interaction with the surface is only through 
the deuterium atoms that point towards the surface. We take the z-axis as the 
surface normal. In this case the surface PES is given by 

A Nd 

^surf—E^, ( 3 ) 
D i=l 

where Njj is the number of deuteriums that points towards the surface, a— 1.0726 
atomic units and A=6A127 Hartree. These parameters are chosen to give the 
same repulsion as the PES that has been used in an MCTDH wavepacket 
simulation of CH 4 dissociation. [26] 

If we write V sur { in terms of normal mode coordinates, then we obtain for one 
deuterium pointing towards the surface 

V surf = Ae- aiXl e~ a2X2 e- a ' iXz e- aiX \ (4) 

where A as above, and a's as given in Table 2. X 2 , X 3 and X A correspond all 
to a\ modes of the C 3 „ symmetry (see Fig. 1). There is no coupling between 
the modes X 5 to X w in the V snr f part of the PES, which are all e modes of 
the symmetry. 

For two deuteriums we obtain 
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Y f — j\_ e ~ a l X l g-"2^2 g-O^Xj g-«5^5 



(5) 



g/3 3 X 7 - /3 3 X 8 g/% Xg - /3 5 Xio 



with A again as above, a's and /3's as given in Table 2. X 2 , X 3 , X A and X 5 
correspond all to a,\ modes of C 2v - Xj, X$, X$ and Xiq correspond to b\ and 
b 2 modes of C 2v . X 6 corresponds to the a 2 mode of C 2v and has no coupling 
with the other modes in V^ ur f. 

For three deuteriums we obtain 



V sulf = A e -^X le -a 2 X 2e -a 3 X 3e -a 4 X 4 



(6) 



X ■ 



g/3i X 5 g/3 2 X 6 — ff s X 7 g/3 4 Xg g/3 5 Xg - j3% X 10 



_|_ g/3iX5g-/32^6g— 03X7 g-^Xg ^f3r,Xg gfleXio 



+ e 



-2/3iX 5 2/33X7-2/35X9 



with A again as above, a's and /3's as given in Table 2. X2, X3 and X4 
corresponds to a\ modes in the symmetry (see Fig. 1). Because these last 
six coordinates correspond to degenerate e modes of the C 3 „ symmetry, the (3 
parameters are not unique. 



2.1.2 An anharmonic intramolecular potential 

Even though we do not try to describe the dissociation of methane in this 
and our previous paper, we do want to determine which internal vibration 
might be important for this dissociation. The PES should at least allow the 
molecule to partially distort as when dissociating. The harmonic PES does 
not do this. A number of changes have therefor been made. The first is that 
we have describe the C-D bond by a Morse PES. 

A 

1 2 



V Mo rs C = D e J2[l-e^ Ar '} , (7) 



where D e = 0.1828 Hartree (the dissociation energy of methane in the gas- 
phase) and Ari the change in bond length from the equilibrium distance. 
7 was calculated by equating the second derivatives along one bond of the 
harmonic and the Morse PES. If we transform Eq. (7) back into normal mode 
coordinates, we obtain 

4 2 

^Morse = D e ^ [l - e 7 "* 2 e 7 * 3 * 3 e 7 * 4 * 4 e 7 * 7 * 7 e 7 * 8 * 8 e 7 * 9 * 9 e 7 *- 10 * 10 ] , (8) 
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with D e as above. 7's are given in Tables 3 and 4. Note that, although we have 
only changed the PES of the bond lengths, the 1S4 umbrella modes are also 
affected. This is because these modes are not only bending, but also contain 
some changes of bond length. 

The new intramolecular PES now becomes 

^intra ^harm ^Morse ^corr; (9) 



where Vh arm is as given in Eq. (2) and V corr is the quadratic part of VMorse, 
which is already in Vharm- 



2.1.3 Intramolecular potential with weakening C-D bonds 

When the methane molecule approach the surface the overlap of substrate 
orbitals and anti-bonding orbitals of the molecule weakens the C-D bonds. 
We want to include this effect for the C-D bonds of the deuteriums pointing 
towards the surface. We have redefined the VM orse given in Eq. (8) and also 
replace it in Eq. (9). A sigmoidal function is used to switch from the gas 
phase C-D bond to a bond close to the surface. We have used the following, 
somewhat arbitrary, approximations, (i) The point of inflection should be at 
a reasonable distance from the surface. It is set to the turnaround point for 
a rigid methane molecule with translation energy 93.2 kJ/mol plus twice the 
fall-off distance of the interaction with the surface, (ii) The depth of the PES 
of the C-D bond is 480 kJ/mol in the gas phase, but only 93.2 kJ/mol near 
the surface. The value 93.2 kJ/mol corresponds to the height of the activation 
barrier used in our dissociation. [26] (iii) The exponential factor is the same as 
for the interaction with the surface. 



If we transform to normal-mode coordinates for the particular orientations, 
we then obtain 



v we , k = D e J2w t 



1=1 



■y e 7i2^2g7i3^3g7i4^"4g7i7^7g7i8^8g7i9^9g7i,10^"l0 



1 2 



\ (10) 



where W{ = 1 for non-interacting bonds and 



1 _L Qp-aiXi+u) 



for the interacting bonds pointing towards the surface. a± is as given in Table 
2, 7's are given in Tables 3 and 4, Q = 1.942 • 10" 1 and co = 7.197. 
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2.1.4 Intramolecular potential with elongation of the C-D bonds 



A weakened bond generally has not only a reduced bond strength, but also an 
increased bond length. We include the effect of the elongation of the C-D bond 
length of the deuteriums that point towards the surface due to interactions 
with the surface. We have redefined the VM orse given in Eq. (8) and also replace 
it in Eq. (9) for this type of PES. We have used the following approximations: 
(i) The transition state, as determined by Ref. [30] and [64], has a C-H bond 
that is 0.54 A longer than normal. This elongation should occur at the turn 
around point for a rigid methane molecule with a translation energy of 93.2 
kJ/mol. (ii) The exponential factor is again the same as for the interaction 
with the surface. 

If we transform to normal-mode coordinates for the particular orientations, 
then we obtain 

4 2 

Vrfuft = D e J2 [l - e 7l2X2 e 7 ' 3X3 e 7l4X4 e 7 ' 7X7 e 7l8X8 e 7 ' 9X9 e 7 '- loXl0 exp^e - " 1 * 1 ]] ,(12) 
i=i 



where a,\ is as given in Table 2, 7's are given in Tables 3 and 4. For orientation 
with one deuterium towards the surface we obtain; S\ = 2.942 • 10 2 and £2 = 
S3 = S4 = 0, with two deuteriums; Si = S 2 = and £3 = S4 = 1.698 ■ 10 2 , 
and with three deuteriums; Si — and S2 — S3 — £4 — 2.942 • 10 2 . 



2.2 Initial States 



The exact wave-function of a D-dimensional system, is expressed in the MCTDH 
approximation by the form 

*mctdh(9i, • ••,?£>;*) = J2 c ni...n D (t)ipS(qi]t) . . . ip^ (q D ; t) . (13) 

n\...n D 



All initial states in the simulations start with the vibrational ground state. 
The initial translational part ip^ is represented by a Gaussian wave-packet, 



^\Xi) = (27ra 2 )- i/4 exp 



2\-l/4 



(Xi - x y< 

Aa 2 



+ tPiXi 



(14) 



where a is the width of the wave-packet (we used a = 320.248 atomic units), 
Xq is the initial position (we used Xo = llcr, which is far enough from the 
surface to observe no repulsion) and Pi is the initial momentum. Since we 
used mass-weighted coordinates the Gaussian wavepacket are identical for CD 4 
and CH 4 . We performed simulations in the energy range of 32 to 96 kJ/mol. 



7 



We here present only the results of 96 kJ/mol (equivalent to P\ = —0.2704 
atomic units), because they showed the most obvious excitation probabilities 
for VMorse- We used seven natural single-particle states, 512 grid points and a 
grid-length of 15a for the translational coordinate. With this grid-width we 
can perform simulation with a translational energy up to 144 kJ/mol. 

Gauss- Hermite discrete- variable representations (DVR)[65] were used to rep- 
resent the wavepackets of the vibrational modes. We used for all simulations of 
CD 4 the same number of DVR points as for CH 4 , which was 5 DVR points for 
the i>2 bending modes and 8 DVR points for the v± umbrella, u 3 asymmetrical 
stretch, and v\ symmetrical stretch mode for an numerical exact integration, 
except for the simulations with V^hift, where we used 16 DVR points for the v\ 
symmetrical stretch mode, because of the change in the equilibrium position. 

Also the same configurational basis was used for both isotopes. We did the 
simulation with one bond pointing towards the surface in eight dimensions, 
because the v 2 bending modes X 5 and X 6 do not couple with the other modes. 
We needed four natural single-particle states for modes X 2 , X 3 and X 4 , and 
just one for the others. So the number of configurations was 7 1 • 4 3 • l 4 = 448. 
The simulation with two bonds pointing towards the surface was performed in 
nine dimensions. One of the v 2 bending mode [X§) does not couple with the 
other modes, but for the other mode X 5 we needed four natural single-particle 
states. The number of configurations was 7 1 -4 4 - 1 4 = 1792, because we needed 
the same number of natural single-particle states as mentioned above for the 
other modes. We needed ten dimensions to perform the simulation with three 
bonds pointing towards the surface. We used here one natural single-particle 
state for the modes X 5 to X w and four natural single-particle states for X 2 
to X 4 , which gave us 7 1 • 4 3 • l 6 = 448 configurations. 



3 Results and Discussion 

3.1 Excitation probabilities and structure deformation of CD 4 

The scattering probabilities for CD4 are predominantly elastic, as we also 
found in our previous simulations of CH4 scattering. [38] The elastic scattering 
probability is larger than 0.99 for all orientation of the PESs with V^rse and 
Kvcak at a translational energy of 96 kJ/mol. For the PES with V^hift we observe 
an elastic scattering probability of 0.981 for the orientation with one, 0.955 
with two and 0.892 with three deuteriums pointing towards the surface. This is 
lower than we have found for CH4, which is 0.956 for the orientation with three 
hydrogens pointing towards the surface and larger than 0.99 for the others. 
The higher inelastic scattering probabilities of CD 4 was expected, because the 
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force constants ki of CD 4 are decreased up to 50% with respect to those of 
CH 4 and the translational surface repulsion fall-off differs only little. 

When we look at the excitation probabilities at the surface for the PES with 
Vyiorse and V wea k, then we observe generally an increase for CD 4 compared with 
CH 4 , except for the z/ 4 umbrella mode in the orientation with two bond point- 
ing towards the surface. Relevant differences in the structure deformations are 
observed only in the bond angles, which are increased for CD 4 in the orienta- 
tions with one and three bonds pointing towards the surface. The bond angle 
deformation of the angle between the bonds pointing towards the surface in 
the orientation with two bonds pointing towards the surface is decreased for 
CD 4 . We observe again that the PES with V wcak gives larger structure defor- 
mations than the PES with \4iorse, but the differences are smaller for CD 4 
than CH 4 . 

For the PES with V^hift we do not observe this effect on the bond angle defor- 
mation. The bond angle deformation for the orientation with two and three 
deuteriums pointing towards the surface is the same as for CH 4 and it is just 
0.1° less for the bond angle at the surface side in the orientation with one deu- 
terium pointing towards the surface. The excitation probabilities (see Table 5) 
for the vi bending and z/ 4 umbrella modes become higher for all orientations 
for CD 4 , which is necessary for getting the same bond angle deformations as 
CH 4 . 

The changes in the bond distances for the orientations with one and two 
bonds pointing towards the surface is for CD 4 almost the same as for CH 4 . 
For the orientation with three bonds pointing towards the surface, we found 
that the maximum bond lengthening of the bonds on the surface side was 
0.032 A less for CD 4 than CH 4 . We also found that the bond shortening of the 
bond pointing away from the surface is 0.010 A more for CD 4 . These are only 
minimal differences, which also only suggest that the bond deformation for 
CD 4 has been influenced slightly more by the z/ 3 asymmetrical stretch mode 
than the V\ symmetrical stretch mode. The observed excitation probabilities 
for these modes do not contradict this, but are not reliable enough for hard 
conclusions because of their high magnitude. It is also not clear, beside of this 
problem, what they really represent. Is the excitation caused by a different 
equilibrium position of the PES at the surface in a mode or is it caused by 
extra energy in this mode? To answer these questions we decided to do an 
energy distribution analysis during the scattering for both isotopes. 
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3.2 Energy distribution in CH4 and CD4 



The energy distribution analysis is performed by calculating the expectation 
values of the important term of the Hamiltonian H for the wave- function ty(t) 
at a certain time t during the scattering of CD 4 and CH 4 for all presented 
orientations in this and our previous paper. [38] We will present here only the 
results of the PES with V^hift, because it is the only model PES for which the 
energy distribution analysis is relevant for the discussion of the dissociation 
hypotheses later on. 

We can obtain good information about the energy distribution per mode by 
looking at the kinetic energy expectation values (^(t)\Tj\^(t)) per mode j 
(see Table 6), because the kinetic energy operators Tj have no cross terms like 
the PESs have. When we discuss the kinetic energy of a mode we normally 
refer to the a\ mode of the 03^, or C 2v symmetry, because in these modes we 
have observed the highest excitation probabilities and the change in kinetic 
energy in the other modes is generally small. 

By looking at the expectation values of some terms of the PES (t)\V term \ty (t)) 
(see Table 7), we obtain information about how the kinetics of the scattering 
is driven by the PES. The V^ ur f PES [see Eqs. (4), (5) and (6)] is the surface 
hydrogen/deuterium repulsion for a given orientation. Vharm^) and Vh arm (i/ 4 ) 
[see Eq. (2)] are the pure harmonic terms of the intramolecular PES of the 
ai modes in the C 3v and C 2v symmetry corresponding to a u 2 bending and z/ 4 
umbrella modes, respectively. The pure harmonic correction terms of V corr [see 
Eq. (9)] are included in them. Vb on d(-R U p) and Vbond(-Rdown) are the potential 
energy in a single C-H or C-D bond pointing respectively towards and away 
from the surface, and they give the expectation value of one bond term of V^hift 
[see Eq. (12)]. All given expectation values are the maximum deviation of the 
initial values, which effectively means the values at the moment the molecule 
hits the surface. 

The largest changes in expectation values are, of course, in the kinetic energy 
of the translational mode. The translational kinetic energy does not become 
zero as we should expect in classical dynamics. The loss of translational kinetic 
energy is primary absorbed by the V sur f terms of the PESs. The expectation 
values of the V sur { terms show the ability of the hydrogens or deuteriums to 
come close to the metal surface, since in real space their exponential fall-offs 
are the same for both isotopes. For a rigid molecule the sum of the translational 
kinetic energy and V^ ur f should be constant, so all deviations of this sum have 
to be found back in the intramolecular kinetic energy and other PES terms. 

We observe that both the minimum in the translational kinetic energy and 
the maximum in the V^ ur f terms were higher for CH 4 than CD 4 , so we have to 
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find more increase in energy in the intramolecular modes and PES terms for 
CD 4 than CH 4 . We indeed do so and that can be one of the reasons we found 
higher inelastic scatter probabilities for CD 4 for the PES with V^hift • 

For the orientations with one and two bonds pointing towards the surface we 
observe a large increase of the kinetic energy in the v$ asymmetrical stretch 
mode. If we compare this with the excitation probabilities, we find that the 
kinetic energy analysis gives indeed a different view on the dynamics. For the 
orientation with two bond pointing towards the surface we have found for both 
isotopes very high excitation probabilities in the v x and i/ 3 stretch modes. We 
know now from the kinetic energy distribution that for the v\ symmetrical 
stretch mode the high excitation probability is caused by the change of the 
equilibrium position of the v\ mode in the PES and that for the v% stretch 
mode probably the PES also has become narrower. 

For the orientation with three bonds pointing towards the surface we also 
obtain an large increase of the kinetic energy of the v?, asymmetrical stretch 
mode, but we also find an even larger increase in the kinetic energy of the 
V\ symmetrical stretch mode. The total kinetic energy was extremely large, 
because the kinetic energy of the translational mode becomes also much larger 
than for the other orientations. Because of this the V^ ur f terms had to be around 
twice as low as for the other orientations. 

All Vbond(-Rup) terms become lower compared to the initial value, especially in 
the orientation with two bond pointing towards the surface. In the orientation 
with one bond pointing towards the surface, the Vbond(-Rdown) term became 
higher. This is caused by the repulsion of V^ ur f in the direction of the bond. 
The increase of this PES term value is higher for CD 4 than CH 4 . 

In the orientation with three bond pointing towards the surface we also observe 
a higher Vbond(-Rdown) value, with also the highest increase for CD 4 . In relation 
with the somewhat shorter bond distance for the i?down of CD 4 compared with 
CH 4 , which was also a bit lower compared with the other orientations, we 
know now that the hydrogens and especially the deuterium have problems in 
following the minimum energy path of the PES with V^hift during the scattering 
dynamics. This leads to higher kinetic energy in the vibrational modes, which 
results in more inelastic scattering. 

The Vharm(^2) term increases in respect to the initial value, but not as much 
as the increase of the Vh arm (^ 4 ) is for the orientation with two bonds pointing 
towards the surface. The values of Vharm^) for CD 4 are even higher than for 
CH 4 . We observe also a larger increase of the kinetic energy in the z/ 4 umbrella 
mode for CD 4 than for CH 4 . So although there is somewhat more energy 
transfer to the vibrational modes for CD 4 than CH 4 , this extra vibrational 
energy is absorbed especially in the z/ 4 umbrella mode of CD 4 . 
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3.3 Dissociation hypotheses 



We like now to discuss some possible implications of the scattering simulation 
for the isotope effect on the dissociation of methane. In our previous paper we 
have already drawn some conclusions about the possible reaction mechanism 
and which potential type would be necessary for dissociation. [38] We found the 
direct breaking of a single C-H bond in the initial collision more reasonable 
than the splats model with single bond breaking after an intermediary Ni-C 
bond formation as suggested by Ref . [4] , because the bond angle deformations 
seems to small to allow a Ni-C to form. From the simulations with CD4 
we can draw the same conclusions. The PES with Khift gives the same angle 
deformations for both isotopes, which is not sufficient for the splats model. The 
other potentials give higher bond angle deformations for the orientation with 
three deuteriums pointing towards the surface. If the Ni-C bond formation 
would go along this reaction path, then the dissociation of CD 4 should be even 
more preferable than CH4, which is not the case. So we only have to discuss 
the implication of the scattering simulation for the dissociation probabilities 
of CH 4 and CD 4 for a direct breaking of a single bond reaction mechanism. 
This reaction mechanism can be influenced by what we will call a direct or an 
indirect effect. 

A direct effects is the expected changes in the dissociation probability between 
CH 4 and CD 4 for a given orientation. Since we expect that we need for disso- 
ciation a PES with an elongation of the bonds pointing towards the surface, 
we only have to look at the isotope effect in the simulation for the PES with 
V^hift for different orientations to discuss some direct effect. It is clear from our 
simulations that the bond lengthening of CD 4 is smaller than CH 4 for the ori- 
entation with three bonds pointing towards the surface. If this orientation has 
a high contribution to the dissociation of methane, then this can be the rea- 
son of the higher dissociation probability of CH 4 . In this case our simulations 
also explain why Ref. [27] did not observe a high enough isotope effect in the 
dissociation probability of their simulation with CH 4 and CD 4 modelled by a 
diatomic, because we do not observe a change in bond lengthening between 
the isotopes for the orientation with one bond pointing towards the surface. 

The orientation with three bonds pointing towards the surface is also the 
orientation with the highest increase of the total vibrational kinetic energy 
for the PES with V^hift, because the energy distribution analysis shows besides 
an high increase of the kinetic energy in the v 3 asymmetrical stretch mode 
also an high increase in the v x symmetrical stretch mode. Since vibrational 
kinetic energy can be used effectively to overcome the dissociation barrier, the 
orientation with three bonds indicates to be a more preferable orientation for 
dissociation. Moreover the relative difference in kinetic energy between both 
isotopes is for the v x stretch mode larger than for the z/ 3 stretch mode. If the 
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kinetic energy in the v\ stretch mode contributes significantly to overcoming 
the dissociation barrier, then it is another explanation for the low isotope 
effect in Ref.[27]. 

An indirect effect is the expected changes in the dissociation probability be- 
tween CH4 and CD4 through changes in the orientations distribution caused 
by the isotope effect in the vibrational modes. This can be the case if the 
favourable orientation for dissociation is not near the orientation with three 
bonds pointing towards the surface, but more in a region where one or two 
bonds pointing towards the surface. These orientations do not show a large 
difference in deformation for the PES with V^hift- We can not draw immediate 
conclusion about the indirect effect from our simulations, since we did not 
include rotational motion, but our simulation show that an indirect isotope 
effect can exist. For the PES with "l4iorse in the orientation with three bonds 
pointing towards the surface, we observe that CD 4 is able to come closer to the 
surface than CH 4 . So this rotational orientation should be more preferable for 
CD 4 than for CH 4 . On the other hand, if the PES is in this orientation more 
like Kshift the dissociation probability in other orientation can be decreased for 
CD 4 through higher probability in inelastic scattering channels. 

So for both effects the behaviour of the orientation with three bonds point- 
ing towards the surface seems to be essential for a reasonable description of 
the dissociation mechanism of methane. A wavepacket simulation of methane 
scattering including one or more rotational degrees of freedom and the vibra- 
tional stretch modes will be a good starting model to study the direct and 
indirect effects, since most of the kinetic energy changes are observed in the 
stretch modes and so the bending and umbrella modes are only relevant with 
accurate PESs. Eventually dissociation paths can be introduced in the PES 
along one or more bonds. 

Beside of our descriptions of the possible isotopes effect for the dissociation 
extracted of the scatter simulations we have to keep in mind that also a tunnel- 
ing mechanism can be highly responsible for the higher observed isotope effect 
in the experiment and that a different dissociation barrier in the simulations 
can enhance this effect of tunneling. 



4 Conclusions 

The scattering is in all cases predominantly elastic. However, the observed in- 
elastic scattering is higher for CD 4 compared with previous simulation on CH 4 
for the PES with an elongated equilibrium bond length close to the surface. 
When the molecule hits the surface, we observe in general a higher vibra- 
tional excitation for CD 4 than CH 4 . The PES with an elongated equilibrium 
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bond length close to the surface gives for both isotopes almost the same defor- 
mations, although we observe a somewhat smaller bond lengthening for CD 4 
in the orientation with three bonds pointing towards the surface. The other 
model PESs show differences in the bond angle deformations and in the distri- 
bution of the excitation probabilities of CD 4 and CH 4 , especially for the PES 
with only an anharmonic intramolecular potential. 

Energy distribution analysis contributes new information on the scattering 
dynamics. A higher transfer of translational energy towards vibrational kinetic 
energy at the surface results in higher inelastic scattering. The highest increase 
of vibrational kinetic energy is found in the asymmetrical stretch modes for 
all orientations and also in the v\ symmetrical stretch mode for the orientation 
with three bonds pointing towards the surface, when the PES has an elongated 
equilibrium bond length close to the surface. 

Our simulations give an indication that the isotope effect in the methane 
dissociation is caused mostly by the difference in the scattering behaviour of 
the molecule in the orientation with three bonds pointing towards the surface. 
At least multiple vibrational stretch modes should be included for a reasonable 
description of isotope effect in the methane dissociation reaction. 
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Fig. 1. The a\ vibrational normal modes in the C^, v symmetry; v\ symmetrical 
stretch (X2), ^3 asymmetrical stretch (X4), and z/4 umbrella (X3). 
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Table 1 

Overview of the relations between the mass- weighted coordinates Xf, the force con- 
stants ki (in atomic units) for CD4, the designation, and the symmetry in T^, Cs v 
and Civ- 



i 


ki 






designation 


T d 


C% v 


C 2 '0 


1 








translation 


t 2 


01 


ai 


2 


8.897 


10" 


-5 


v\\ symmetrical stretch 


ai 


a\ 




3 


2.008 


10" 


-5 


1/4; umbrella 


t 2 


ai 


ai 


4 


1.060 


10- 


-4 


1/3; asymmetrical stretch 


t 2 


ai 


ai 


5 


2.447 


10" 


-5 


v 2 ; bending 


e 


e 


m 


6 


2.447 


10" 


-5 


u 2 ; bending 


e 


e 


a 2 


7 


2.008 


10" 


-5 


1/4; umbrella 


h 


e 


h 


8 


2.008 


10" 


-5 


1/4; umbrella 


t 2 


e 


b 2 


9 


1.060 


10" 


-4 


v%\ asymmetrical stretch 


t 2 


e 


h 


10 


1.060 


10' 


-4 


1/3; asymmetrical stretch 


t 2 


e 


b 2 



Table 2 

a and (3 values (in atomic units) of V SUT f for CD4 with one, two or three deuteriums 
pointing towards the surface [see Eqs. (4), (5) and (6)]. 



one two three 





5.617 


10" 


-3 


5.617 


10" 


-3 


5.617 


10" 


-3 


a 2 


8.882 


10" 


-3 


5.128 


10" 


-3 


2.960 


10" 


-3 


a 3 


4.703 


10" 


-3 


-4.614 


10" 


-3 


-7.720 


10" 


-3 


04 


-1.353 


10" 


-2 


-5.103 


10- 


-3 


-2.295 


10" 


-3 


a 5 








-7.252 


10" 


-3 








Pi 














4.187 


10" 


-3 
















7.252 


10- 


-3 










4.659 


10" 


-3 


2.196 


10' 


-3 
















3.804 


10" 


-3 










4.212 


10" 


-3 


2.295 


10- 


-3 


fa 














3.439 


10" 


-3 
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Table 3 

7 values (in atomic units) of VMorse f° r CD4 with one and three deuteriums pointing 
towards the surface [see Eq. (8)]. 



one 




thx66 




value 


712,722,732,742 




712, 722, 732, 742 


7 69Q 


• 1 n -3 


713, -3723, -3733, 


-3743 


-713, 3723 , 3733 , 3743 


1.397 


■ icr 3 


714, -3724, -3734, 


-3744 


-714,3724,3734,3744 


-1.454 


• 1(T 2 


717, 718, 719, 71,10, 728, 72,10 


717, 718, 719, 71,10, 728, 72,10 




0.0 


727, -2737, -2747 




-727, 2737, 2747 


1.318 


• 10~ 3 


738, —748 




738, ~748 


-1.114 


• 10~ 3 


729, -2739, -2749 




-729, 2739, 2749 


-1.371 


• 10- 2 


73,10, -74,10 




-73,10,74,10 


1.187 


• lO- 2 



Table 4 

7 values (in atomic units) of VMorse for CD4 with two deuteriums pointing towards 
the surface [see Eq. (8)]. 

two value 

712, 722, 732, 742 7.629 • 10~ 3 

713, 723, -733, -743, 717, -727, 737, -747, 718, "728, —738, 748 -8.070 • 10 -4 
714, 724, -734, -744, 719, -729, 739, -749, 7i,io, -72,10, -73,io, 74,io 8.396 • 10~ 3 



Table 5 

Excitation probabilities at the surface, at an initial translational energy of 96 kJ/mol 
and all initial vibrational states in the ground state, for the intramolecular PES 
with elongation of the C-D bonds [see Eq. (12)] in the a\ modes of the C3,, and C2 V 
symmetry, with one, two or three deuteriums pointing towards the surface. These 
modes are a ui(a±) symmetrical stretch, a v<i{e) bending, a ^(£2) asymmetrical 
stretch, and a ^4(^2) umbrella. In parenthesis is the irreducible representation in T^ 
symmetry. 



orientation 


v\{a\) stretch 


V2 (e) bending 


^3 (£2) stretch 


^4(^2) umbrella 


one 


0.460 




0.910 


0.174 


two 


0.792 


0.092 


0.830 


0.495 


three 


0.868 




0.756 


0.387 
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Table 6 

Expectation values of the kinetic energy per mode in mHartree for CH4 and CD4, 
at an initial translational energy of 96 kJ/mol and all initial vibrational states in 
the ground state, for the intramolecular PES with elongation of the C-H/D bonds 
[see Eq. (12)] in the a\ modes of the C% v and C^u symmetry, with one, two or three 
deuteriums pointing towards the surface. These modes are a vi{a\) symmetrical 
stretch, a v<i{e) bending , a ^3^2) asymmetrical stretch, and a ^4^2) umbrella. In 
parenthesis is the irreducible representation in T^ symmetry. 



isotope 


orientation 


translation 


v\(a\) stretch 


v-2, (e) bending 


^3(^2) stretch 


^4 (£2) umbrella 


CH 4 


initial 


36.57 


3.30 


1.75 


3.39 


1.50 




one 


16.76 


3.56 




4.53 


1.51 




two 


14.59 


3.50 


1.79 


4.67 


1.57 




three 


20.53 


5.32 




4.39 


1.58 


CD 4 


initial 


36.57 


2.33 


1.24 


2.52 


1.12 




one 


16.17 


2.61 




4.09 


1.18 




two 


14.00 


2.78 


1.27 


4.05 


1.27 




three 


20.06 


4.37 




3.80 


1.28 



Table 7 

Expectation values of the potential energy terms in mHartree for CH4 and CD4, at 
an initial translational energy of 96 kJ/mol and all initial vibrational states in the 
ground state, for the intramolecular PES with elongation of the C-H/D bonds [see 
Eq. (12)]. V sw { is the total surface hydrogen repulsion; T4axm(^2) and 14 ar m(^4) ar e 
the harmonic terms of the intramolecular PES of the a\ modes in the C^y and C2 V 



symmetry corresponding to a v^ie) bending and ^(£2) umbrella modes respectively 
in the T^ symmetry. Vb on d(-Rup) an d Vbond(-^down) are the potential energy in a 
single C-H/D bond pointing respectively towards and away from the surface. 


isotope 


orientation 


^surf 


Vharm^) 


Vharm(l / 4) 


Vbond(-Rup) 


^bond (-Kdown) 


CH 4 


initial 


0.00 


1.75 


1.50 


3.39 


3.39 




one 


18.20 




1.87 


3.25 


3.85 




two 


18.55 


2.18 


4.01 


2.75 


3.45 




three 


9.22 




2.94 


3.00 


3.74 


CD 4 


initial 


0.00 


1.24 


1.12 


2.48 


2.48 




one 


17.94 




1.89 


2.43 


3.44 




two 


18.45 


1.68 


4.52 


2.29 


2.74 




three 


8.71 




3.49 


2.28 


3.21 
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